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Measurements of the elliptic flow, V2, of identified hadrons (it , K , K®, p, p, <j>, A, A, H~, 

E + , Cl~ , in Au+Au collisions at ^/snn = 7.7, if. 5, 19.6, 27, 39 and 62.4 GeV are presented. 
The measurements were done at mid-rapidity using the Time Projection Chamber and the Time- 
of-Flight detectors of the STAR experiment during the Beam Energy Scan program at RHIC. A 
significant difference in the V2 values for particles and the corresponding anti-particles was observed 
at all transverse momenta for the first time. The difference increases with decreasing center-of- 
mass energy, ^snn (or increasing baryon chemical potential, hb) and is larger for the baryons as 
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compared to the mesons. This implies that particles and anti-particles are no longer consistent with 
the universal number-of-constituent quark (NCQ) scaling of V2 that was observed at ^/snn = 200 
GeV. However, for the group of particles NCQ scaling at (hit- mo)/n q > 0.4 GeV/c 2 is not violated 
within ±10%. The V2 values for <f> mesons at 7.7 and 11.5 GeV are approximately two standard 
deviations from the trend defined by the other hadrons at the highest measured pt values. 

PACS numbers: 25.75.Ld, 25.75.Nq 



I. INTRODUCTION 

One of the main goals of the heavy ion collision pro- 
gram at the Relativistic Heavy Ion Collider (RHIC) fa- 
cility is to produce a state of deconfined quarks and 
gluons, called the Quark-Gluon Plasma (QGP), and to 
study its properties. An experimental way to understand 
the formation of the QGP is by varying collision ener- 
gies and studying observables as a function of collision 
centrality, transverse momentum, pr, and rapidity, y. 
This also allows one to study the structure of the QCD 
phase diagram. With these goals, the Beam Energy Scan 
(BES) program was started in the years 2010 and 2011 
at RHIC P| where Au+Au collisions were recorded at 
= 7.7, 11.5, 19.6, 27, 39 and 62.4 GeV. This pa- 
per reports the azimuthal anisotropy of identified par- 
ticles produced in collisions at BES energies, measured 
using the STAR detector at RHIC. 

The azimuthal anisotropy of produced particles is one 
of the most widely studied observables. In non-central 
heavy ion collisions, the overlap region of the colliding 
nuclei is almond-shaped and perpendicular to the plane 
defined by the impact parameter vector and the beam 
axis. This plane is called the reaction plane. Due to fi- 
nite number fluctuations of participating nucleons in re- 
actions with the same centrality, the geometric symmetry 
plane in each event is not necessarily the same as the re- 
action plane, and is often called the participant plane. 
This plane is defined by the nucleons which participated 
in the reaction Q. In a hydrodynamic approach with 
local thermalization, the initial spatial anisotropy and 
subsequent interactions among the constituents result in 
pressure gradients that are larger in the direction of the 
participant plane compared to out of this plane. This 
results in an azimuthal anisotropy of the momenta of the 
produced particles Q. The second harmonic parameter, 
t>2, of the Fourier decomposition of the azimuthal particle 
distribution relative to the event plane is called the ellip- 
tic flow [1, Q . It is experimentally measured using final- 
state particle momenta. The event plane is an approxi- 
mation to the participant plane. The elliptical anisotropy 
with respect to the event plane is not necessarily equal to 
the elliptic flow with respect to the participant plane. In 
the literature, the magnitude of the second flow harmonic 
is called V2 whether this quantity is calculated from the 
participant (or reaction) plane or the event plane. The 
expansion of the system and subsequent decrease of the 
spatial anisotropy leads to a self-quenching process for 
V2, thereby making it a sensitive probe of the early stage 
of heavy ion collisions 0, 0] • 



Using the data from the top RHIC energy of 200 GeV 
several interesting observations related to V2 have been 
reported in the past decade pj-[l2|. Large values of 
the elliptic flow were found to be compatible with ideal 
hydrodynamic calculations [l3|-[l5j or viscous hydrody- 
namic calculations [ljl [ijj with a small shear viscosity 
to entropy density ratio. At low transverse momentum 
(pr < 2 GeV/ c), a mass ordering of the v 2 values was ob- 
served [1S| - [20| ] , which could readily be understood within 
a hydrodynamic framework. At intermediate px values 
(2 < pt < 6 GeV/c), a Number-of-Constituent Quark 
(NCQ) scaling of V2 for the identified hadrons was ob- 
served. This observation, coupled with the compara- 
ble values of the elliptic flow measured for multi-strange 
hadrons (0 and 5) and light quark hadrons, was used 
to conclude that the relevant degrees of freedom in the 
systems formed at the top RHIC energy are quarks and 
gluons [3.l20l-l23l|. It was also concluded that a substantial 
amount of v 2 was generated during the partonic stage of 
these collisions. This was further corroborated by com- 
paring the measurements to model calculations with and 
without partonic interactions. 

It is generally expected that the system will spend less 
time in the partonic phase as the beam energy is low- 
ered, and that at the lowest BES energies the system 
might not reach the QGP regime. In such a scenario, it 
is expected that NCQ scaling of V2 of produced parti- 
cles would be broken [24|]. Furthermore, with decreasing 
beam energy, the baryon chemical potential of the sys- 
tem at chemical freeze-out increases. These aspects could 
lead to new trends in the identified hadron V2 in the BES 
program at RHIC which was performed at the BES ener- 
gies with unmatched statistics and particle identification 
capabilities. In this paper, the STAR measurements of 
the beam energy and px dependence of the elliptic flow, 
V2, at mid-rapidity for 7r ± , K ± , K®, p, p, <f), A, A, H~, 
5 , 0~ and il in minimum bias (0-80%) Au+Au col- 
lisions are presented. The corresponding results for the 
inclusive charged particles were reported in Ref. [25| . 

This paper is organized as follows. Section |H] gives 
a brief description of the experimental setup and the 
event and centrality selection. In Sections IIIIHIV1 andfVl 
the various particle identification methods, the event 
plane reconstruction, v 2 signal extraction, and system- 
atic uncertainty estimation are discussed. In Section IVI1 
the energy-and-momentum-dependent V2 results are pre- 
sented. In Section IVII1 comparisons to models are 
discussed. Finally, the summary is presented in Sec- 
tion Enu 
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II. EXPERIMENTAL SETUP 

The Solcnoidal Tracker At RHIC (STAR) is a multi- 
purpose experiment at the RHIC facility at Brookhaven 
National Laboratory. It consists of a solenoidal magnet 
and an array of detectors for triggering, particle identifi- 
cation, and event categorization. A detailed description 
can be found in Ref. [26[. The primary detectors used 
for the present results are summarized in the following 
subsections. 



A. Time-projection chamber (TPC) 

The TPC has a full azimuthal, (j), coverage and a pscu- 
dorapidity, r/, acceptance of —1.8 < r\ < 1.8 [27}. The 
TPC is split into two halves along the beam direction by 
a central membrane. A maximum of 45 hit points per 
track can be reconstructed within the TPC radius lim- 
its of 0.5 < r < 2 m. The primary collision vertex of 
an event is fitted using the reconstructed particle tracks. 
For ~1000 such tracks, a primary vertex resolution of 350 
lim can be achieved. The primary vertex position is used 
in a subsequent track refitting for particles like 7T, K and 
p to improve the momentum resolution. The relative mo- 
mentum resolution for pions is ^3% at Pt = 1 GeV/c. 
The specific energy loss (dE/dx) information, also pro- 
vided by the TPC, can be used for particle identification 
(c/. Section EDI. 



B. Time-of-Flight (TOF) 

The time-of-flight system is based on Multi-gap Re- 
sistive Plate Chambers (MRPCs) and was fully installed 
in STAR in the year 2010 [Hj|. The system has an in- 
trinsic timing resolution of ~85 ps. It covers the full 
azimuth and a pseudorapidity range of —0.9 < r/ < 0.9. 
The matching efficiency of a TPC-reconstructed track 
to an MRPC cell is ^90%, which results in a total ef- 
ficiency (acceptance x efficiency) of ^65%. The particle 
mass-squared, m 2 , can be calculated using the measured 
time-of-flight and the reconstructed momentum from the 
TPC. Examples of the m 2 distributions are shown in Sec- 
tion [nil 



C. Trigger and event selection 

In the years 2010 and 2011, Au+Au collisions at the 
six energies, ^snn, of 7.7, 11.5, 19.5, 27, 39 and 62.4 
GeV were measured. The minimum bias trigger condi- 
tion for all six energies was based on a coincidence of the 
signals from the Zero Degree Calorimeters (ZDC), Vertex 
Position Detectors (VPD), and/or Beam-Beam Counters 
(BBC). Most of the triggered events at the lowest beam 
energies did not originate from Au+Au collisions, but 



were rather Au plus beam pipe (or other material) col- 
lisions. This was the result of the large beam emittance 
at the lowest beam energies. The radius of the beam 
pipe is 3.95 cm. The background due to these "fixed tar- 
get" events was efficiently removed in the present anal- 
ysis by requiring that the primary vertex position was 
within a radius r of less than 2 cm. The z-position of the 
primary vertices was limited to the values listed in Ta- 
ble. HI These values depend on the offline ^-vertex trigger 
conditions which were different for the different energies. 
These vertex cuts were studied and optimized during the 
data-taking using the online vertex reconstruction per- 
formed by the High Level Trigger (HLT). 

To remove pile-up events, it was required that at least 
two tracks from the primary vertex were matched to the 
cells of the time-of-flight detector. Furthermore, an ex- 
tensive quality assurance of the events was performed 
based on the mean transverse momenta, the mean ver- 
tex position, the mean interaction rate, and the mean 
multiplicity in the detector. Run periods were removed 
if one of those quantities was several a away from the 
global mean value. The accepted number of minimum 
bias events for each of the six energies are also listed in 
Table. |U 



D. Centrality definition 

The centrality selection of the events was chosen to 
be 0-80% of the total reaction cross section. The cen- 
trality definition was based on the uncorrected multiplic- 
ity distribution, dN cvts /dN^ w , of reconstructed charged 
particle tracks within a pseudorapidity range of \r)\ <0.5. 
The distributions for all energies can be accurately de- 
scribed by a 2 component model calculation Q as shown 
in Fig. [T] Some of the most peripheral events were not 
recorded due to trigger inefficiencies. This results in a sig- 
nificant difference between the measured dN evts /dN*™ 
and the Glauber Monte Carlo (MC) simulation. To cor- 
rect for this effect, the ratio of the simulation to the data 
was used as a weighting factor for the particle yields. 
The correction drops from a maximum of 30% to 5% by 
the 70% most central bin, and is negligible for the more 
central data. In addition to the trigger inefficiency cor- 
rections, two additional corrections were also applied to 
account for the z-vertex dependent inefficiencies. These 
corrections treated the acceptance and detector ineffi- 
ciencies and the time dependent changes in dN 0Vts / dN^ w 
resulting, e.g., from minor changes in the trigger config- 
uration. 



III. PARTICLE IDENTIFICATION AND 
SIGNAL EXTRACTION 

Particle identification in the STAR experiment can be 
done in multiple ways. Long-lived charged particles, e.g. 
7T, K and p, were directly identified and reconstructed, 
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FIG. 1. (Color online) The uncorrected multiplicity, , distribution of reconstructed charged particles per unit pseudo- 
rapidity interval at mid-rapidity for the six different beam energies. The filled black points depict the measured data and a 
Glauber Monte Carlo simulation is overlayed as the solid curve. Three different centrality classes are indicated by the different 
shaded regions. 



TABLE I. The total number of minimum-bias (MB) events used, and the z-vertex acceptance, for the different energies. 



'sjvat (GeV) MB events (10 b ) z-vcrtcx range (cm) 

TT~ O - [-70,70] 

11.5 11.7 [-50,50] 

19.6 35.8 [-70,70] 
27 70.4 [-70,70] 
39 130.4 [-40,40] 

62.4 67\3 [-40,40] 



within a pseudorapidity range of \rj\ < 1, using the time- 
of- flight information and/or the specific energy loss in 
the TPC depending on the reconstructed track momen- 
tum. For weakly-decaying particles, e.g. A and S, the 
invariant mass technique and topological reconstruction 
methods were used. They are reconstructed within a 
rapidity range of \y\ < 1. The cleanest event-by-event 
particle identification is obtained at the lowest momenta 
and/or when using tight topology cuts. Statistical signal 
extraction methods were used to obtain the yields of the 
particles at higher momenta. Up to momenta of ~1.5 
GeV/c a clean separation of 7T, K and p was obtained 
when combining the information from the TPC and TOF 
detectors. At higher momenta, the it and K signals begin 
to overlap. The protons still can be separated event-by- 



event up to ~3.0 GeV/c by using the time-of-flight infor- 
mation alone. Figure [2] shows the mean specific energy 
loss, (dE/dx), in the TPC and the mass-squared from 
the TOF as a function of the momentum. The proton, 
pion, and kaon dE/dx bands merge for momenta above 
~1 GeV/c. The separation in m 2 of tt, K and p at a 
beam energy of ^/snn = 19.6 GeV is shown for three 
different momentum ranges in Fig. [3J 

In order to avoid fake tracks in the TPC and to im- 
prove the average momentum and energy loss resolution, 
the following track quality cuts were applied: the number 
of total hit points was larger than 15, and the ratio of the 
number of reconstructed hits to the maximum possible 
number of hits for each track was larger than 0.52. The 
momentum of each particle was limited to 0.15 < p < 10 
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FIG. 2. (Color online) The mean specific energy loss, 
(dE/dx), of reconstructed tracks within a pseudorapidity 
range of \n\ < 1 in the TPC a), and the mass-squared, m 2 , as 
a function of momentum b) . The Bichsel functions [29[ used 
to determine the no- pa rticie values (cf. Eq. JT}) are shown in 
a) as the dashed curves. The horizontal dashed lines in b) 
correspond to the nominal particle masses of it, K and p. 



GeV/c. The deviation in units of <7 par ti c ie of (dE/dx) of 
a particle species from its theoretical energy loss, calcu- 
lated with a Bichsel function [2{|, can be expressed as, 



^-Cparticle 



(X In 



dE\ 

~dx / 



I 



particle 



(-) 
\dx/ Bi 



Bichsel 



(i) 



The distribution of er p articie is nearly Gaussian for a given 
momentum and is properly calibrated to be centered at 
zero for each particle species with a width of unity. 



A. Signal extraction for ir , K , p, and p 

Protons and anti-protons are identified primarily using 
the TOF mass-squared information. To suppress rem- 
nant contributions from pions and kaons, an additional 
dE/dx cut of \no p \ < 3 was applied. At low transverse 
momenta (pr <2 GeV/c), the separation of protons rela- 
tive to pions and kaons was sufficient such that all protons 
in a range of ~ 3<r around the center of the na p distribu- 
tion are counted. At high pr, the tails on the low mass 
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from 0-80% central Au+Au collisions at the beam energy of 
19.6 GeV. Three different momentum ranges are shown. 



range of the proton distributions were excluded to avoid 
contamination from pions and kaons. Thus, the m 2 cuts 
increased with the transverse momentum, pt- 

For the analysis of tt^ and K^, a new technique was 
employed to extract the yields for each pt bin. This was 
based on a transformation of the combined TOF m 2 and 
TPC dE/dx nov information. The goal of this trans- 
formation was to have a maximal separation between 
kaons and pions by transforming to a new set of vari- 
ables x,y(ncr w ,m 2 ) such that the widths of the parti- 
cle peaks in x and y were identical and for which the 
pion and kaon peaks were aligned with the horizontal 
axis. Each particle was described by two two-dimensional 
(2x2D) Gaussians (a;,y(ncr 7r , m 2 )), where the first Gaus- 
sian fits the peak and the second Gaussian shares the 
same position as the first, but the width was larger to 
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a) 27 GeV, 0-80%, 2.2<p r <2.4 GeV/c 




FIG. 4. a) The mass-squared, m 2 , versus nov and b) x^yina^, m 2 ) (see Eqs. (J2]| to (|6])) distributions for 2.2 < pr < 2.4 GeV/c 
from 0-80% central Au+Au collisions at 27 GeV. The black dashed contour lines in b) depict the result of a simultaneous fit 
with three 2x2D Gaussians. The diagonal dashed line depicts a cut to remove the remaining proton contamination (see text), 
c) The projected distribution to the x(na^, m 2 ) axis. The red solid curve shows the projection of the 2x2D Gaussian fits, d) 
The same as c), but after the 2x2D Gaussian of the protons was removed. The red solid line shows the sum of the two ID 
Gaussian fits. The fit range is indicated by the two vertical dashed lines. 



account for the broad tail. The tt, K and p peaks of 
the m 2 vs. na^ distributions are fit simultaneously, in- 
dividually for each px bin with three 2x2D Gaussians. 
The non-Gaussian tails of the tt, K and p peaks along 
both axis were excluded from the fits. The resulting fit 
parameters, widths a(m 2 )(Tr) and a(na„)(7r), and peak 
positions, {j,(iti 2 )(tt, K) and /i(nov)(7T, K), were used to 
first normalize the m 2 axis to the na axis and then to 
perform a transformation which consists of a shift and a 
rotation. The transformations are listed in Eqs. ([2] to [5]), 



/scale = a^ina^ja^m 2 ), 



(2) 



a 



X/ 



V {fi K (nan ) - {na^ ) ) / / sca i c J 
(na^ - /i w (n<J w )) // sca ic, (4) 



m 



x(n<Jis, m 2 ) 
y{na^, m 2 ) 



cos(a) — sin(a) 
sin(a) cos(a) 



(5) 



(6) 



Figure [4] a) shows an example of the m 2 versus na n 
distributions and frame b) shows the new x,y(ncr n , m 2 ) 
distribution after the transformation for an intermediate 
transverse momentum range of 2.2 < px < 2.4 GeV/c. 
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The protons were treated as background in the 7r ± , 
analysis and were removed first. For this, the distribu- 
tions in the new x,y{na v , m 2 ) frame were fitted with 
three 2x2D Gaussians in a way analogous to that de- 
scribed above. To stabilize the procedure, several itera- 
tions were performed. In the final fit, only the area 2.5 
u(x, y) away from the pion and proton peak positions was 
considered. The fit range and the resulting fit arc shown 
as a dashed box and dash-dotted contour lines, respec- 
tively, in Fig.Ub). The projection to the x(na v , m 2 ) axis 
of the data and the fit are shown in Fig. 2] c) . The data 
can be well described for all px bins with the fit function, 
which allows one to subtract the 2D proton fit function 
from the distribution. In addition to the fit subtraction, 
a mass-squared cut of m? < 0.65 (GeV/c 2 ) 2 was applied, 
as shown in Fig.@]b) as a diagonal line. The latter cut re- 
moves the remnant non-Gaussian tails from the protons. 
The corresponding result after the proton subtraction is 
shown in Fig. 3] d) . This distribution was fitted with two 
Gaussians (2xlD) to extract the pion and kaon yields. 
The goal of this representation, the increased separation 
power between the pions and kaons along the transformed 
horizontal axis, was reached. 

B. The signal extraction for <f>, A, A, K°, H~, H + , 
f2~, and Q + 

Short-lived weakly-decaying particles, generically 
called V° particles, such as A, <fi and 3, decay into a pair 
of oppositely charged particles and were reconstructed 
using the invariant mass technique. The combinatorial 
background from uncorrclated particles was reduced by 
a direct identification of the daughter particles using the 
specific energy loss and/or mass-squared (m 2 ) informa- 
tion and selection criteria based on the topology of the 
specific decay. Depending on the particle species and the 
magnitude of the background, no cuts of ±2 or ±3 were 
applied to the normalized dE/dx of the daughter par- 
ticle tracks. Since the time-of-flight information is only 
available for about 65% of the tracks within the accepted 
pscudorapidity range of— 1 < < 1, a general cut on 
the mass-squared m 2 , as for dE/dx, was not applied. In- 
stead, a cut on m 2 was only applied if the time-of-flight 
information for the track was available or the misidenti- 
fication rate at a certain momentum range, when using 
only the dE/dx information, was large. For most of the 
V° particles, the combinatorial background can be effi- 
ciently reduced with topology cuts as will be described 
below. In these cases, a lack of TOF information was 
compensated for by using tighter topology cuts. For the 
<j) meson, the time-of-flight information was always re- 
quired for daughter tracks at higher momenta, typically 
at p > 0.65 GeV/c where the dE/dx information alone 
was insufficient to remove the bulk of the misidentified 
tracks. In general, a 3er cut on the m 2 distributions of 
the particles was applied. At higher momenta, the 7T, K 
andp distributions begin to overlap. Here, tighter and/or 



asymmetric cuts were used. 

For the topological reconstruction of V° particles, ge- 
ometrical information on the decays was also used, e.g. 
the primary and secondary/tertiary decay vertex po- 
sitions, the distance of closest approach (dca) of the 
daughter particles to the primary vertex, the dca of the 
mother particlc(s) to the primary vertex, and the dca 
between the daughter tracks. This information was de- 
termined from the helix paramctcrizations of the TPC 
reconstructed tracks. The following topology cuts were 
applied: 

• dca between daughter tracks (primary and sec- 
ondary daughters in case of 3 and fi), 

• dca between daughter tracks and primary vertex, 

• dca between mother particle and primary vertex, 

• dca between A candidate and primary vertex (for 
3 and f2), and 

• distance between primary and secondary (tertiary 
in case of 3 and Q) vertex. 

A cut on the invariant mass of 1.108 < M(p, ir) < 1.122 
GeV/c 2 was applied to enhance the A and A candidates 

for the 5~, 3 + , and Q~ , fl + analyses. The particle 
identification and topology cuts were systematically opti- 
mized for the best significance by varying several tens of 
thousands of cut combinations for each particle species. 

The misidentification of the daughter particles, which 
is more probable at the higher momenta, can result in 
an additional correlated background. Such a correlated 
background, for example from the A, can appear in the 
7r + 7r~ {K®) invariant mass distribution if the proton was 
misidentified as a n + . Such a correlated background 
does not create a peak in the invariant mass distribu- 
tion of the particles of interest since the daughter par- 
ticle masses are chosen to be the nominal ones {e.g. tt 
mass instead of proton mass), but it appears as a broad 
distribution which can significantly affect the signal ex- 
traction. To remove this correlated background, addi- 
tional invariant mass spectra with identical track combi- 
nations, but different daughter mass values, e.g. (p, tv~ ) 
and (7r + ,7r~) were investigated. The background was 
removed by applying invariant mass cuts on the corre- 
sponding unwanted peaks in the misidentified invariant 
mass distributions. Usually, the correlated background 
from particle misidentification increases with the pr val- 
ues of the mother particle. 

The remaining uncorrclated combinatorial background 
was described and later subtracted with the mixed event 
technique. Event classes were defined to mix only events 
with similar global properties; the classes consisted of 9 
centrality ranges, 14 z-vertex ranges, and 10 event plane 
angle ranges. The event buffer depth varied between 3 
and 15. The mixed event distributions were normalized 
at least 3<r away from the mass peak on both sides. The 
mixed event distributions so obtained were in excellent 
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FIG. 5. Examples of the invariant mass distributions at ^/sjvjv = 62.4 GeV for 0, A'", A, A, H", S , Q~ , and f2 . The 
combinatorial background is described by the mixed event technique which is shown as a grey shaded histogram. 



agreement for all particle species and energies with the 
combinatorial background shown in Fig. H] for 0-80% cen- 
trality Au+Au collisions at v / s]v"/v=62.4 GeV for 0, A, A, 

K®, S~, S + , f2~, and f2 + . The correlated background 

to the left of the E _ , S + , 0~ and il + peaks in Fig. [5] 
is a result of a self-correlation between the three daugh- 
ter particles. In the case of the S, two ir mesons with 
the same charge are in the final state and both com- 
binations (A(p,7ri)+7T2 and A(p,7T2)+7ri) result in similar 
invariant mass values. These wrong combinations appear 
as a bump structure to the left of the true peak. The two 
structures were separated by an invariant mass cut. In 
the case of the f2, a double misidcntification of the ir and 
K resulted in a similar effect. 



IV. EVENT PLANE RECONSTRUCTION 

The event plane (EP) is obtained from the angles of 
the reconstructed particles and the beam line. It is an 
estimate of the participant plane which is defined by the 
participating nucleons in the collision. The event plane 
was reconstructed using the flow of the measured parti- 
cles, as discussed in Ref. @. To achieve the best reso- 
lution for t>2, the second harmonic event plane angle 'J 2 
was calculated as: 



\& 2 = tan~ 



J2i w i sin(2fa) 
J2i ^ cos(20j) 



(7) 



where fa is the azimuthal angle of particle i and Wi is 
its weight. The weight in units of GeV/c was chosen to 
be linear with px up to 2 GeV / c and then constant at a 
value of 2 for higher momenta. Only those particles with 



a momentum between 0.15 and 5 GeV/c, \rj\ < 1, dca 
< 1 cm and having more than 15 hits in the TPC were 
used for this calculation. Two different event planes were 
reconstructed: one using all of the reconstructed tracks 
in the TPC ("full TPC" method) and one using only 
those tracks in the opposite pseudorapidity hemisphere 
to the particle track of interest ("77-sub" method). In 
the full TPC case, self-correlations were avoided by re- 
moving the particle of interest from the tracks used for 
the event plane reconstruction. In the 77-sub method, 
an additional pseudorapidity gap of ±0.05 was applied 
to reject some tracks for the event plane reconstruction. 
In general, the 77-sub method reduces the effect of "non- 
flow," which includes the decay of resonances to several 
charged daughterjparticles, Hanbury-Brown Twiss corre- 
lations, and jets [8j. However, the resolution is lower and 
therefore the correction to obtain u 2 is larger. 

An azimuthally non-homogeneous acceptance or effi- 
ciency of the detectors can introduce a bias in the event 
plane reconstruction which would yield a non-uniform 
^2 angle distribution in the laboratory coordinate sys- 
tem. To flatten the \&2 distribution, the recentering or 
fa weight methods, in combination with the shift method, 
were used In the fa weight method, a track- by-track 
correction is applied. It is based on the </> angle distribu- 
tions, dN / d<j){z,pT, 77, t, q) which were determined for five 
z-vertex ranges, four pt ranges, six 77 ranges, and for both 
charge signs, q. Furthermore, the distributions were de- 
termined for different real time, t, periods during the data 
collection, each of which spanned approximately one day. 
Each track used for the event plane reconstruction was 
weighted in the 0-wcight method with the inverse value 
of the corresponding value of the dN/ dfaz, px, 77, t, q) dis- 
tribution. Large gaps in the 4> angle distribution cannot 
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FIG. 6. The event plane resolution for the full TPC event plane (circles), and the 77-sub event plane (stars), as a function of 
the centrality for two different and independent flattening methods. 



be corrected with this method. The ^-weight corrected 
event plane angles will be denoted as ^2,<p- 

The rcccntcring method applies a correction on an 
event-by-event basis and is therefore more robust in the 
case of acceptance holes. The numerator and denomina- 
tor of Eq. ([7]) can be used to define the vector 

n — 1 \^ I w i cos(20 4 ) \ , , 

V™- N 2^ I tu iB iii(2&) J ' W 

where N is the number of tracks used for the event plane 
reconstruction in each event. In order to get a uniform 
\&2 angle distribution, this Q-vector must be centered at 
(0,0). To achieve this objective, the average of the Q- 
vector over many events was subtracted event- by- event: 

Qrc Qraw ^Qraw^ ■ (9) 

These averaged Q-vectors were determined for ten z- 
vertex ranges and for each pseudorapidity hemisphere 
in a real-time dependent manner, and were then 
parametrized as a function of the event multiplicity. The 
new ^2,rc angles were then calculated from the corrected 
Q- vectors. 

If the ^2 angle distribution was not flat after the <fi- 
wcight or the recentering corrections, an additional cor- 
rection with the shift method was used to force the ^2 
angle distribution to be flat Q. A shift angle ^2, shift was 



calculated event- by- event for each event plane method in 
the following way: 

*2,shift = -c 2 cos(2* 2 ) + s 2 sin(2* 2 ) 

+ 0.5(-c 4 cos(4* 2 ) + s 4 sin(4* 2 ))- (10) 

The C2,4,S2,4 parameters were obtained from fits to the 
averaged (^-weight or recentering corrected ^2 angle dis- 
tributions. The shift-corrected event plane angle v I , 2,corr 
was obtained as: 

*2,corr = *2,rc,0 + #2,shift, (H) 

where ^2,rc,<f> is the recentering corrected event plane an- 
gle. After the shift correction a flat Vf^.corr distribu- 
tion for all energies and event plane methods was thus 
achieved. 

To calculate the event plane (EP) resolution, indepen- 
dent sub-samples of randomly selected tracks (full TPC) 
or tracks in independent pseudorapidity hemispheres (77- 
sub) were used [8[ . Figure [5] shows the event plane res- 
olution for the four different reconstructed event plane 
types and the six beam energies. The event plane res- 
olution is used below (cf. Section IV Bp to correct the 
observed ?j§ bs signals. The event plane resolution is ap- 
proximately proportional to the flow coefficient times the 
square- root of the multiplicity It decreases with de- 
creasing beam energy due to the lower particle multi- 
plicities. It has a maximum for each beam energy at 
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about 30% centrality. For more peripheral events, the 
relatively low multiplicity is responsible for the decreas- 
ing resolution whereas for more central events the small 
flow signal is responsible. The 0- weight corrected EP has 
a slightly smaller resolution compared to the recentering 
method which could be connected to the smaller number 
of centrality bins used for the 0-weight correction. At 
62.4 GeV, a significantly larger difference between the 
two correction methods is observed compared to all of 
the other beam energies. This is due to a missing TPC 
sector during the collection of the 62.4 GeV data. The 
resulting gap in the ^2 angle distribution cannot be fully 
corrected with the </>-weight method as described above. 

In general, the 77-sub method has a smaller EP resolu- 
tion compared to the full TPC method. This is mainly 
due to the factor of ^2 fewer tracks used for the EP re- 
construction in the former. For the most central collisions 
and the lowest energies, 7.7 and 11.5 GeV, a similar EP 
resolution for the two methods is observed. This might 
be an indication of a strong negative non-flow signal at 
the lower energies when the full TPC method is used. 
The negative non-flow, which originates primarily from 
resonance decays, results in an anticorrelation between 
the random sub-events used for the EP resolution calcu- 
lation. The 77-sub EP method reduces the non-flow by 
using spatially independent regions in the TPC. There- 
fore, in the following only the results based on the 77-sub 
EP method will be presented. 



V. v 2 SIGNAL EXTRACTION AND 
SYSTEMATIC UNCERTAINTIES 

The azimuthal emission pattern of the particles relative 
to the event plane can be decomposed into a Fourier sum 
of cosine terms: 



oc 1 + 2 y v n cos [n( 



d(</>-* 



n>l 



(12) 



where 4> is the azimuthal angle of the particle, ^> m is the 
event plane angle, v n is the Fourier coefficient of har- 
monic n, and m is the harmonic of the event plane 
In the following, only the elliptic flow coefficient v%, will 
be considered, 
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FIG. 7. (Color online) Two examples of the V2 signal extrac- 
tion for (f> mesons at 39 GeV in the transverse momentum 
range of 0.8 < pr < 1.0 GeV/c. The event plane method a) 
and the invariant mass method b) give almost identical re- 
sults, a) The <f> — ^2 data points are reflected at n/2. A fit 
with Eq. (| 13p to the data obtained by integrating the fit is 
shown as a solid black line. The red dashed line shows the fit 
result to the data obtained by counting the particles in each 
bin. b) The solid black curve is the fit from Eq. (|15p . The 
dashed red curve is the signal part of that equation and the 
dashed blue curve is the background part. 



ocl + 2u 2 cos[2(0-tf 2 )]- (13) 



A. Event plane and invariant mass methods 

Two techniques were used to calculate 7J 2 : the event 
plane method and the invariant mass method 0. The 
latter method was used in addition only for the <p mesons. 
Both methods give identical results, but are technically 
different, especially when the signal-to-background ratio 



is small. In the standard event plane method, the par- 
ticles were first identified, then their yields were deter- 
mined as a function of the relative angle <j) — \&2- In the 
invariant mass method, the mean values (cos [2(0 — \I/ 2 )]) 
were calculated as a function of invariant mass and then 
the correlation at the invariant mass peak of interest was 
isolated after the background subtraction. 

For the event plane method, the 1*2 coefficients were 
obtained by fits to the yield distributions with Eq. (|13l) . 
An example of such a fit is shown in Fig. [7] a). For most of 
the particles, the yields were determined in two different 
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ways: by counting the particles in bins within reason- 
able mass ranges and / or by integrating a fit to the corre- 
sponding mass distribution. The fits can have the form 
of a Gaussian or a Breit-Wigner distribution. A Brcit- 
Wigner distribution was used for the <f> mesons, and a 
Gaussian was used for the A and particles. The yields 
were only determined in the <p — ^ 2 range of to tt/2. 
In Fig. [7] a), the reflected data points are also shown for 
reference. 

For the invariant mass method, the mean values (...) 
(see also Ref. [30j ) were by definition the v 2 values of 
the analyzed particles (c/. Eq. (TT3)) ). Since the back- 
ground cannot be distinguished from the signal on an 
event- by- event basis, the resulting v 2 value was the sum 
of signal and background as formulated in Eq. ([M]) . The 
f 2 IS+Bs (M inv ) can be decomposed into a signal and a 
background term as shown in Eq. (TT5]). Each term is 
multiplied by a statistical weight which was extracted 
from the same event invariant mass and the combinato- 
rial background distributions. The background elliptic 
flow u^ s (Minv) was parameterized with the polynomial 
defined in Eq. (|16|) . Figure [7] b) shows an example of 
an invariant mass fit with the total fit (black solid line), 
the signal term (red dashed line), and the background 
term (blue dashed line). In this particular case, the v 2 e 
is nearly identical to u 2 ls which results in a monotonic 
distribution around the signal region. 

v S 2 is+Bs (M iav ) = (cos [2(0 - * 2 )] Mfnv ) (14) 

V2 Sig+Bs (Mnv) = vfZ-^i—(M inv ) 
Sig + Bg 

+ v^(M inv ) 68 (A/ inv ) (15) 
Sig + Bg 

w Bg (Mnv) = po + PlM inv + p 2 M? nv 

+ PsML (16) 

The invariant mass method was tested for various parti- 
cle species and directly compared to the results from the 
event plane method. For particles with large signal-to- 
background ratios in the invariant mass distribution, for 
instance A and H, no systematic differences were found. 
The present results are generally based on the event plane 
method. However, both methods were evaluated only for 
the 4> meson, which shows a significantly lower signal-to- 
background ratio compared to all other particles. Small 
differences between the two methods were taken into ac- 
count in the systematic uncertainties. 

B. Event plane resolution correction for 0—80% 

The event plane resolution was calculated for nine cen- 
trality bins as shown in Fig. [6] For the integrated 0-80% 
ccntrality bin, a new method was used to correct the 
observed v% bs si gnals. The yields of the reconstructed 
particles were weighted event-by-event with the inverse 



event plane resolution for the corresponding ccntrality 
bin. This ensured a correction which was not biased by 
the bin width. The v 2 signals were normalized with the 
mean inverse event plane resolution for the 0-80% cen- 
trality bin. A detailed description of the method can be 
found in Ref. [30l |. 



C. Systematic uncertainties 

The systematic uncertainties were evaluated by vary- 
ing the methods and parameters used to determine the 
event plane angles and particle yields. For the V° parti- 
cle analyses, e.g A and H, twenty different combinations 
of the topology cuts listed in Section fill Bl were applied. 
Each of these topology cut combinations had a signifi- 
cance similar to the reference cuts which were optimized 
for the best significance. The same number of combi- 
nations were used for the <f> meson analysis, but in this 
case it was the K ± particle identification cuts that were 
varied, such as the nax range. 

In addition to the variation of particle identification 
cuts, two methods to extract the v 2 values and two ways 
to determine the particle yields were used as described in 
Section [V Al The it and K analyses depend primarily on 
the initial fit parameters and fit ranges as pointed out in 
Section Till Al The parameters were varied in combination 
with two values of the proton separation m? cuts, result- 
ing in a total of eighteen different combinations. In the 
case of the proton analysis, three different combinations 
of dca and na p cuts were studied. 

The point-by-point systematic uncertainties, which 
consist of the variations of the particle identifications cuts 
and the two methods of signal extraction, were evaluated 
for all combinations (40 for all V° particles and the <j> me- 
son, 36 for 7r and K, and 6 for protons) by calculating 
the root- mean-squared value for each data point. For all 
energies, both flattening methods for the event plane an- 
gle (which were described in Section HV]) were compared. 
The mean value of the point-by-point differences between 
the two methods was defined as the global systematic un- 
certainty for each particle species. The mean point-by- 
point systematic uncertainties varied for p, tt and K in 
a range of 0.0001-0.001, and for V° particles and the <j> 
meson in a range of 0.0005-0.007. The mean global sys- 
tematic uncertainties for all particles were in the range 
of 0.0005-0.003. In addition to these studies of the sys- 
tematic uncertainties, independent analyses for most of 
the particle species were performed. A cross check to the 
previously published 62.4 GcV data (22|, where slightly 
different methods were used, indicated an excellent agree- 
ment to the present results within the statistical errors. 

The data were not corrected for feed-down contribu- 
tions. By varying the dca cuts for all particle species, 
the feed-down contributions were already partly included 
into the systematic uncertainties described above. Previ- 
ous studies showed [2l| that feed-down is only significant 
for pions below pt=0A GeV/c. For other particles, the 
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resulting feed-down contributions to the v 2 values are 
negligible. 

Non-flow contributions were studied for the six beam 
energies by comparing different methods to extract v 2 for 
inclusive charged hadrons (25j . The four- particle cumu- 
lant i>2{4} strongly suppresses non-flow contributions. It 
has been shown that the difference between v 2 (77-sub) and 
v 2 {4} is about 10-20% for 19.2, 27, and 39 GeV and de- 
creases with decreasing beam energy. In the following we 
did not treat such non-flow contributions as systematic 
errors. 



VI. RESULTS 

The i>2 results corrected for the event plane resolution 
in 0-80% central Au+Au collisions are presented. All 
results are based on the 77-sub event plane method de- 
scribed above. The x-axis values of the data points are 
always placed at the p^-wcightcd mean values within the 
bin limits. The statistical errors are indicated as straight 
vertical lines, the point-by-point systematic uncertain- 
ties are indicated either as shaded bands or with square 
brackets, and the global systematic uncertainties are in- 
dicated as a horizontal shaded band on the horizontal 
axis. For plots with several v 2 distributions, only the 
statistical errors are shown. 



Elliptic flow as a function of transverse 
momentum 



Figure [5] shows the energy dependence in v 2 (pt) for 
, K^, p, p, A, A, <fi, Kg, S~ and S + . A similar trend 
of v 2 (pt) for all particles is observed. The v 2 increases 
with pt up to 1.5 GeV / c and reaches a maximum value at 
higher beam energies of about 0.15 for mesons and ~0.2 
for baryons within the measured px range. The max- 
imum values decrease with decreasing energy to about 
0.07 for kaons and pions and ~0.1 for protons at ^s^n 
= 7.7 GeV. It should be noted that the v 2 (pt) decreases 
for higher p T values in ^/snn = 200 GeV Au+Au col- 
lisions (3l| . The negative anti-proton v 2 at low pt and 
at v / saF/v = 11.5 GeV could be due to absorption in the 
medium (32|. A more precise picture of the energy de- 
pendence can be obtained from v 2 (px) ratios. In order 
to define a reference, the 39 GeV data points were fitted 
with the following equation: 



fv 2 {n) 



an 



1 + e^^VT/n-b)/c 



(17) 



where a, b, c and d are fit parameters and n is the 
constituent-quark number of the particle [33|. The cor- 
responding ratios are shown in Fig. [9] A non-trivial pt 
dependence is observed from the ratios of the v 2 (pr) val- 
ues to the fits. The ratios are close to unity and nearly 
independent of pt for y/s^N > 19-6 GeV. Below 19.6 
GeV, the ratios decrease (below unity) with increasing 



Pt for ^ , p, A and K + . At these energies the ratios are 
below unity, but are independent of pt for p, A, K® and 
K~. 

As expected, the ^2{??-sub}-energy ratios of charged 
particles [25[ follow the same trends as presented here 
for the identified hadrons. The trends for the more abun- 
dantly produced particles presented here (p, it, K) differ 
from those obtained using the inclusive charged hadron 
four-particle cumulant, t>2{4} [25|. For the lower beam 
energies, the values of the ratios for the inclusive charged 
hadron «2{4} increase with increasing px- As will be 
discussed below, the difference might be due to non-flow 
contributions and flow fluctuations. For ir^ , p, A and 
K + , the ratios for all energies are close to unity at low 
Pt and deviate with increasing px, whereas the ratios for 
p, A, and K~ seem to be independent of pr for all 
energies. 

At low transverse momenta, a mass ordering was ob- 
served it Au+Au collisions at 200 GeV [3l( . Lighter par- 



ticles had larger v 2 values. This behavior can be qualita- 
tively described by ideal hydrodynamics [34j. In Fig. ITUI 
the v 2 (pt) values in the transverse momentum range of 
0.2 < pt < 1.6 GeV/c for various particle species are 
directly compared. For this selection of particles (p, A, 
E~, Q~, 7T+, K + , AT° and 0), the mass ordering is valid 
for all energies. Only the <f> mesons deviate from this 
general trend at the lower energies. Their v 2 (pt) values 
are slightly smaller compared to all of the other hadrons. 
Starting at 39 GeV, every <p meson v 2 {pt) value is smaller 
than the corresponding value for the heavier A. 

The lower the energy, the smaller is the difference be- 
tween the various particles in v 2 (j>t) at pt < 1.5 GeV / c. 
This could be related to a reduction of radial flow as the 
beam energy decreases. However, no narrowing of the 
spread of v 2 (pt) with beam energy is observed for the 
anti-particles. At lower beam energies, the v 2 {pt) values 
for p and A were significantly smaller than the values for 
their partner particles. The possible physics implications 
due to the differences in particle and anti-particle v 2 (j>t) 
will be discussed in more detail in the next sections. 

In Fig.[TTJ each particle v 2 (pt) is directly compared, if 
possible, to that for its anti-particle. For the mesons 
the anti-particle convention from [35| is used. The 
point-by-point systematic uncertainties are displayed as 
the shaded bands which connect the data points. The 
global systematic uncertainties are shown as the error 
bands along the horizontal axis. Shown are the v 2 (pt) 
for TT+(ud), w-(ud), and K+(us), K°((ds - sd)/\/2), 
K~(us). At the higher energies of 27, 39 and 62.4 GeV, 
the charged pion 7r + and n~ v 2 (pt) values show almost 
identical shapes and amplitudes, as expected from parti- 
cles with the same mass and number of quarks. At lower 
energies, an increasing difference between i>2(7r + ) and 
v 2 {~k~) is observed, where v 2 (ir~) is larger than v 2 (ir + ) 
for all pt values. In the lower rows of each panel in 
Fig. [n] the difference in v 2 (pt) between particles and 
anti-particles is shown. The red line shows a horizontal 
line fit to the Av 2 which will be used below (Section I VI D() 
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FIG. 8. (Color online) The elliptic flow, V2, as a function of the transverse momentum, pr, from 0-80% central Au+Au collisions 
for various particle species and energies. Only the statistical error bars are shown. The black dashed line is a fit to the 39 GeV 
data points with Eq. (|17|l . 



to study the energy dependence of the difference. The fit 
range was varied to estimate the systematic uncertainty 
for the fit and to test the assumption of a constant dif- 
ference as a function of pt- 

The fact that v-2(ir~) is larger than V2(ir + ) could be due 
to the Coulomb repulsion of ir + by the mid-rapidity pro- 
tons or to the chiral magnetic effect in finite baryon den- 
sity matter produced in the collisions [36[. The charged 
kaons show an opposite trend compared to the charged 



pions. The V2{pt) values of K + are larger compared to 
K~ . The size of the difference in V2 and the energy de- 
pendence is comparable to that of the pions. The neutral 
K° approximately follow the trends of the V2(pt) values 
of the R -. 

In contrast to the charged pions and kaons, a signifi- 
cant difference in the V2{pt) values between p{uud) and 
p(uud) was already observed at 62.4 GeV, as shown in 
Fig. Q~2]a). The difference in V2 is nearly constant as a 
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FIG. 9. (Color online) The ratio of the elliptic flow, V2(pr), relative to a fit to the 39 GeV V2(pr) data points for 0-80% central 
Au+Au collisions for various particle species and energies. The error bars are statistical only. The fit and the V2(pr) data 
points are shown in Fig. [8] 
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FIG. 10. (Color online) The elliptic flow, V2(pr), in 0-80% central Au+Au collisions for selected particles a) and anti-particles 
b) (see text), plotted only for the transverse momentum range of 0.2 < pr < 1.6 GeV/c to emphasize the mass ordering at low 
Pt- The fit functions correspond to Eq (|17[) . 



function of pt and, as for the pions and kaons, the dif- 
ference increases with decreasing energy. Compared to 
the kaons and pions, the relative difference is at least 
a factor of three larger. The plots in Fig. [T^] b) show 
the corresponding ^(.Pt) for A(uds) and A(uds). The 
shapes and magnitudes of V2(pt) for all energies are al- 
most identical between p and A and the same between 
p and A. Hence, the difference in V2(j>t) between the 
(anti)A particles and the (anti)protons is observed. It 
appears that the exchange of a u-quarks with an s-quark 
has no influence on the difference in u 2 (pt)- 

Figure IT51 a) shows the V2{pt) of S~(dss) and S + (dss) 



and b) shows the V2(pt) of Q (sss) and ft (sss). Within 
the statistical and systematic uncertainties, S - and S + 
are indistinguishable in V2(pr) at 62.4 GeV. At 39 and 
27 GeV, only a slightly larger v 2 (pt) of S~ with respect 

to S + is observed, whereas at 19.6 and 11.5 GeV the 
difference is significant and comparable to that of the 
protons and A. Due to the larger error bars, no significant 
effect is observed for the and f2 + at any energy. 

As mentioned above, the 4>(ss) meson V2(pt) is of par- 
ticular interest. The hadronic cross section of </> mesons 
is much smaller compared to that of other hadrons [37i — 
|40| . This would result in a smaller V2(pt) for a fireball 
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FIG. 11. (Color online) The elliptic flow, V2, of charged pions a) and kaons b) as a function of the transverse momentum, px, 
for 0-80% central Au+Au collisions. The point-by-point systematic uncertainties are shown by the shaded areas, otherwise 
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FIG. 12. (Color online) The elliptic flow, V2, of p, p a) and A, A b) as a function of the transverse momentum, pr, for 0-80% 
central Au+Au collisions. The point-by-point systematic uncertainties are shown by the shaded areas, otherwise they are 
smaller than the symbol size. The global systematic uncertainties are shown as the shaded horizontal bar. The lower row of 
each panel depicts the difference between a particle and corresponding anti-particle i>2 (pt) with a fit with a horizontal line. 
The red shaded area around each fit shows the combined statistical and systematic fit errors. 



19 



a) 



0.2 



3~ 0.1 



M 0.1 
^ 0.05 


0.2 
j~ 0.1 


.0.02 

< o 

-0.02 



b) 



- 7.7 GeV 

■ 

■ • s + 

" -ft- E -S 
■ 


Au+Au, 0-80%- 
rj-^ub EP 

4 
r 


. i i i i i _ 

-11.5 GeV 

. 

- 


-19.6 GeV 

1 

: aS* T - 


: >t' 


t 


' 1 i i i i 1 i i i i 1 i i i i 1 i i i i 1 " 


; i i i i i i i i i i i i i i i i i i i i i ; 


'-- \ 




Li J 

; '- 

r— 

t ; 












- 27 GeV 

o * - 

8 


i i i i i i i i i i i i i i i i i i i i i _ 
- 39 GeV 

§ '& 

r 

9 

: 


- 62.4 GeV fi 

§ 

* J 




;i i i i i i i i i i i i i i i i i i i i i; 












L ^* 1 j 




12340 12340 1234 

P T (GeV/c) 



> 



> 

< 



0.2 


-0.2 
0.2 



0.2 


-0.2 


-0.1 



7.7 GeV 

- o a 
• n 


Au+Au, 0-80% 
-q-sub EP 


L 1 1 1 1 1 1 

11.5 GeV 


i i i i i i i i i i i i _ 

f : 


- 19.6 GeV 

L.iLL.J 






< 

- I i i i i 1 ■ . i 


, i , , , , i , . , , i - 










[_l i . i i I i i i i I i i i i I i i i i I - 
'- -JtmmtS,- '- 








c. . ^~~¥ .: 


- 27 GeV 

'. A 


t 


1 , , , , | . , , , | , , , , | , . , , | _ 

■ 39 GeV 

* 

§ 


- 62.4 GeV 

$ 






;l i i i i 1 i i i i 1 i i i i 1 i i i i 1: 


; 1 i i i i 1 i i i i 1 i i i i 1 i i i i 1 ; 






t ' ': 


- a . 


: t i x_ : 








i- _ 

c .¥ - 


ymmrntf _ 



12340 12340 1234 

P T (GeV/c) 
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FIG. 15. (Color online) The elliptic flow, V2, of p and p as a function of the transverse momentum, pr, for 0-10% central 
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evolution in the hadron gas phase. The results arc shown 
in Fig. [Til At 19.6 to 62.4 GeV, the typical v 2 (pt) shape 
is seen, whereas at 7.7 and 11.5 GeV, the v 2 values at the 
highest measured pr bins are close to zero. Also, there is 
a significant decrease in the energy dependence of v 2 (px) 
(cf. Fig. [5]) at transverse momenta of about 1.5-2 GeV/c. 

In [4l[ v 2 (pt) studies at -^/saT/v = 17.3 GeV for pions, 
kaons, and strange particles are presented for mid-ccntral 
Pb+Au collisions. Due to the different ccntrality selec- 
tion, a direct comparison was not performed. 



three of these narrower centrality ranges, a significant 
difference between the p and p v 2 (px) values is seen. For 
the most peripheral ccntrality bin (40-80%), the elliptic 
flow is the largest, but the absolute difference Av 2 (pt) is 
smaller compared to the mid-central bin (10-40%) and is 
comparable to the most central bin (0-10%). It is con- 
cluded that Av 2 (pt) shows a clear ccntrality dependence 
for protons and anti-protons, and that the difference in 
v 2 (pt) remains when restricted to narrower centrality 
ranges. 



1. Centrality dependence of proton and anti-proton V2 



B. Elliptic flow as a function of transverse mass 



The elliptic flow shows a strong ccntrality dependence 
which is driven by the changing initial spatial eccen- 
tricity. The present results are an average over a wide 
(0-80%) ccntrality range. Even if the v 2 (pr) values for 
protons and anti-protons would be identical for all col- 
lision centralities, one would observe a difference in the 
v 2 values if the centrality dependency of the production 
rates would be very different. To study this possibility, 
Fig. \U)\ Fig. [T51 and Fig. [T7] show the proton and anti- 
proton v 2 (pt) values for the centrality ranges of 0-10%, 
10-40%, and 40-80% for the six beam energies. In all 



The v 2 values as a function of the reduced transverse 
mass, tot ~ shows a clear splitting between baryons 
and mesons for larger — too values at y / s]v"/v =200 
GeV (T^j. The particle mass, charge, and strangeness 
content are not the driving factors. Only the number 
of constituent quarks separates the results into the two 
branches. This observation is an indication that the re- 
sults are sensitive to the particle internal degrees of free- 
dom, i.e. the quarks in the QGP phase of the collision. 
After hadronization, the flow of the quarks is carried by 
the measured particles. In a coalescence picture, this will 
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FIG. 17. (Color online) Elliptic flow, vi, of p and pas a function of transverse momentum pr for 40-80% centrality Au+Au 
collisions. The point-by-point systematic uncertainties are shown by the shaded areas, while the global systematic uncertainties 
are shown as the shaded horizontal bar. Shown in the lower row of each panel is the difference between a particle and the 
corresponding anti-particle V2(pt) which are fit with a horizontal line. The red shaded area around each fit shows the combined 
statistical and systematic fit error. 



result in the V2 values of the baryons being a factor of 1.5 
larger than the V2 values of the mesons [13[ . Figure [TS] 
shows the ^(tot — mo) values for all six BES energies 
and the same selection of particles a) and correspond- 
ing anti-particles b) as presented above. The baryons 
and mesons are clearly separated in Fig. [T5] a) above 
(triT — mo) > 1 GeV/c 2 . The separation at 7.7 GeV be- 
tween protons and 7r + , K + is significantly smaller than 
that at all of the other energies. The A hypcrons follow 
the meson branch at 7.7 GeV. 



The anti-particles at 39 and 62.4 GeV show a similar 
behavior as the particles, and at all lower energies the 
meson and baryon branches approach each other. At 11.5 
GeV, a difference between the baryons and mesons is no 
longer observed, and at 7.7 GeV the anti-proton and A 
V2(mT—mo) are below the meson branch in the measured 
rriT — mo range. The trend observed is a decrease in the 
baryon-meson splitting in ^(m^ — mo) for (ttt-t — mo) > 
1 GeV/c 2 as the energy is lowered, both for the particle 
and anti-particle groups. 



C. Number-of-constituent quark scaling of V2 

The splitting in vvivriT — mo) between the mesons and 
baryons at transverse mass values above 1 GeV/c 2 shown 
in Fig. [TS] implies a dependence of the V2 values on the 
number of constituent quarks, n q . The NCQ scaling was 
originally predicted for V2{pr) at intermediate transverse 
momenta [42[ ■ A scaling of pr and V2 with l/n q was sug- 
gested. Indeed, the scaled V2 values for all particles at 
200 GeV Au+Au collisions collapse to a common single 
trend at intermediate px values (2fj| - [23l . [43| . This is in- 
terpreted as a possible signature for partonic degrees of 
freedom (quarks and gluons) in the initial stage of the 
system, where most of the elliptic flow develops. This 
scaling should vanish in a hadron gas system at lower 
energies. Thus, the break down of NCQ scaling would be 
a necessary signature for a QCD phase transition from 
partonic to hadronic matter. 

Since particles and anti-particles have the same num- 
ber of quarks, the NCQ scaling transformation of V2 does 
not change their relative separation. This means that the 
difference in V2(j>t) for particles and corresponding anti- 
particles observed in Section IVI Al constitutes a violation 
of this NCQ scaling. Possible physics causes for this dif- 
ference will be discussed below. In the following, NCQ 



23 



a) 



b) 



> 



- 7.7 

; Au+i 

rj-su 


GeV Op *7i + - 
<\u, 0-80% aa nK + " 
b EP Az *K° - 


- 1 1 .5 GeV 

rt. 


. ... i .... i .... i .... i . 

19.6 GeV 










_ I I 1 1 _ 

- 27 GeV 


i i i i I i i i i I i i i i I i i ' i I <['- 
- 39 GeV 


' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I 
62.4 GeVj ^ A : 

jm&T* $ : 









0123 40 123 40 1234 

m T -m (GeV/c 2 ) 




m T -m (GeV/cr) 



FIG. 18. (Color online) The elliptic flow, V2, of 0-80% central Au+Au collisions as a function of the reduced transverse mass, 
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is observed at the higher energies. The splitting becomes smaller at 7.7 GeV. At lower energies, the baryons and mesons are 
consistent with each other within the measured pr range for the particles shown in b). 



scaling will be shown separately for a selection of parti- 
cles and anti-particles. Since a better agreement between 
the different particles (even at low (tot — mo)/n q values) 
is achieved with the (v2 / n q )((rnT — m$)/n q ) scaling com- 
pared to the (v2/n q )(pT/n q ) scaling, Fig. [TO] presents the 
scaled distributions versus (tot — mo)/n q . The corre- 
sponding scaled plots for V2(pr) are shown in Fig. [M] of 
the Appendix. 

The NCQ scaling should only hold in the transverse 
momentum range of 1.5 < pr < 5 GeV/c [3, E3|- For 
the corresponding scaled transverse mass and transverse 
momentum range, a fair agreement for most of the parti- 



cles and energies is observed. Only the <f> mesons deviate 
from the trend at 7.7 and 11.5 GeV, with the maximum 
measured PT/n q value just reaching the lower edge of 
the expected NCQ scaling range. The values deviate 
from those for the other particles and anti-particles at 
the highest (tot — too) values at i/sJjn = 7.7 and 11.5 

GeV by 1.8a and 2.3er, respectively. The values for S + 
at 11.5 GeV were similar and show a relatively small V2 
compared to the other hadrons. This could be related to 
the lower hadronic cross sections of particles containing 
multiple strange quarks. These observations may indi- 
cate that hadronic interactions become more important 
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FIG. 19. (Color online) The Number-of-Constituent Quark (NCQ) scaled elliptic flow, V2/n q versus (bit — mo)/n q , for 0-80% 
central Au+Au collisions for selected particles a) and corresponding anti-particles b). The dashed lines show the results of 
simultaneous fits with Eq. (|17[) to all particles except the pions. 



than partonic effects for the systems formed at collision 
energies < 11.5 GeV [HI, H^]. It is noted that recent 
results of elliptic flow measurements of J/ty mesons at 
y/SNN = 200 GeV also show smaller V2(pt) values com- 
pared to those for other hadrons (47|- However, in this 
representation of NCQ scaling, both particles and anti- 
particles appear to follow the scaling. In the previous 
sub-section, an absence of the baryon-meson splitting of 
vi{m,T — mo) for (rriT — m ) > 1 GeV/c 2 is observed. 

For each energy, simultaneous fits with Eq. (|17[) were 
applied to all particles except the pions, which are biased 
by resonance decays HU . The ratios of the data to the fits 
are shown in Fig. [2D] and for the transverse momentum 



in Fig. [35] of the Appendix. Most of the data points in 
the high transverse momentum range, agree within the 
uncertainties within a ±10% interval around unity. At 
lower values of hit — Too, larger deviations from unity 
are observed. 
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FIG. 20. (Color online) The Number-of-Constituent Quark (NCQ) scaled elliptic flow, vijn q versus (mr — mo)/n q , ratio to the 
fit function (see text) for 0-80% central Au+Au collisions for selected particles a) and corresponding anti-particles b). Most 
of the data points at the larger (mr — mo)/n q values are within a ±10% interval around unity, which is shown as the shaded 
region to guide the eye. Some data points for (f> and S are outside of the plot axis range. 



D. Energy dependence of the particle and 
anti-particle V2 difference 



In this sub-section, the energy dependence of the v 2 
difference between particles X (p, A, S~, ir + , K + ) and 

anti-particles X (p, A, S + , ir~ , K~) is studied. Figurcl2"T1 
a) shows a fit to the Av 2 {pt) values from Figs.HTl[T2l and 
H"3l This difference is denoted in the following as v 2 {X) — 
v 2 (X) and is shown as a function of the beam energy 
\J s nn ■ At 62.4 GeV, the v 2 difference for mesons is close 
to zero, whereas the baryons show a difference of 0.003 to 



0.005. The difference increases for all particle species as 
the energy decreases. It reaches values of about 0.03 for 
A and protons and 0.004-0.005 for kaons and pions at 7.7 
GeV. The baryons show a steeper rise compared to the 
mesons. The pions and kaons show a similar trend, but 
opposite with respect to their charge. Also, the protons 
and A are very similar at all energies. Compared to the 
protons and A, the S show a slightly smaller difference at 
higher energies, but a larger difference at lower energies. 
One should note that the S result at 11.5 GeV covers a 
much smaller px range compared to all of the other data 
points. This could cause additional systematic effects 
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TABLE II. Fit parameters a and b of Eq. (|18[) and the slopes m of the straight line fits shown in Fig. [21] b) for the different 
particle species. The first \ 2 P er Number-of-Degrees of Freedom (NDF) value corresponds to the fit with Eq. (|18[) . the second 



to the straight line fits. Both include statistical and systematic uncertainties. 

Particle a b \ Z /NDF m x T jNDF 

7T+ - tF -0.064 ± 0.025 1.068 ± 0.106 22.5 / 4 -0.0155 ± 0.0005 22.8 / 5 

K+-K- 3.219 ± 10.068 3.104 ± 1.4440 0.4 / 4 0.0018 ± 0.0017 7.3 /5 

p-p_ 0.209 ±0.099 0.9329 ±0.143 1.8 / 4 0.0831 ± 0.0039 2.0 /5 

A — A 0.177 ±0.086 0.896 ±0.139 0.6 / 4 0.0794 ± 0.0040 0.7 /5 

7.363 ± 18.997 2.072 ± 0.825 1.0 / 3 0.0607 ± 0.0210 2.5 /4 



which are not included in the error bars. The difference 
in i>2 (a/sjvjv ) shown in Fig. [21] a) was parametrized with: 

/au 2 Wsnn) = a ■ s^% 2 . (18) 

The fit results of the parameters a and b are listed in 
Table H] 

In Fig. [5T] b) , the i>2 difference is shown as a function 
of the baryonic chemical potential, [1b- A parametriza- 
tion from [48| was used to determine the \ib values for 
each beam energy. Since this parametrization was done 
for the most central collisions, a correction has to be ap- 
plied to take into account the difference to the minimum 
bias collisions. To do this, the measured centrality de- 
pendence of hb for Au+Au collisions at ^snn = 62.4 
and 200 GeV from [49[ was used. The ratio between the 
mean fxs over all centralitics, which is an approximation 
for the 0-80% central fiB values, and the most central 
values is 0.83 ±0.06 for 62.4 GeV and 0.84 ±0.14 for 200 
GeV. The minimum bias [Lb values were calculated for all 
energies by multiplying this factor with the obtained val- 
ues from the parametrization under the assumption that 
these ratios do not change with energy. The resulting [Ib 
values with errors are shown in Fig. 1211 b). 

Each particle data set was fitted with a straight line 
that passes through the origin; the slope parameters, to, 
are listed in Table [IT] A linear increase of the V2 differ- 
ence with [Ib is observed for all particle species from 62.4 
GeV down to 7.7 GeV. Only at 11.5 GeV a 2cr devia- 
tion for the S and at 7.7 GeV a deviation for the kaons 
was found. This linear scaling behavior suggests that 
the baryon chemical potential is directly connected to 
the difference in vi between particles and anti-particles. 

VII. DISCUSSION 

Comparisons of the data to transport and other models 
are described. 



A. Transport model comparisons 

In Fig. [221 the measured elliptic flow of tt^, K^, p and 
p for 0-80% central Au±Au collisions at ^/snn = 7.7, 
11.5 and 39 GeV is compared with model calculations. 
The models used were UrQMD (Ultra-relativistic Quan- 
tum Molecular Dynamics), version 2.3 and AMPT 



(A Multi-Phase Transport), version 1.11 [51j . In order 
to be consistent with the analysis of the data, the num- 
ber of charged particles within the pseudorapidity range 
\r)\ < 0.5 was used for the centrality definition in the 
model calculations. The t>2 values of the (anti-)particlcs 
were calculated relative to the true event plane. The dif- 
ference between the true event (or participant) and the 
reaction plane (as inferred in the data) could bias the re- 
sults as discussed in Ref. [52[ . In total, about one million 
events were generated from each model at each energy. 

The UrQMD model [5(| is based on a microscopic 
transport theory where the phase-space description of the 
collisions plays the central role. It allows for the covari- 
ant propagation of all hadrons on classical trajectories 
including stochastic binary scattering, color string for- 
mation, and resonance decay. This model includes more 
than fifty baryon and forty-five meson species and incor- 
porates baryon-baryon, meson-baryon and meson-meson 
interactions. A comparison of the data with the UrQMD 
model can provide information about the contributions 
to the elliptic flow from the hadronic phase. 

The AMPT model [EH has Glauber-based initial con- 
ditions which are the same as those used in the HIJING 
(Heavy Ion Jet Interaction Generator) (53} model. In this 
model, mini-jet partons are created and scatter before 
they fragment into hadrons. The String Melting version 
(AMPT-SM) of the AMPT model is based on the idea 
that, for energy densities beyond a critical value of about 
1 GcV/fm 3 , strings and partons cannot coexist. There- 
fore, the strings are melted into partons by converting the 
hadrons into their valence quarks. The Zhang's Parton 
Cascade (ZPC) model [5l| was used to describe the scat- 
tering between the quarks. Once the interactions have 
stopped, the partons hadronize through the mechanism 
of parton coalescence. In the default AMPT model, par- 
tons are recombined with their parent string when they 
stop interacting, and the resulting strings are converted 
to hadrons using the Lund string fragmentation model. 
The interactions between the mini-jet partons in the de- 
fault AMPT model and those between partons in the 
AMPT-SM model could give rise to substantial elliptic 
flow. The AMPT-SM calculations would thus indicate 
the contributions to the measured v% from the partonic 
interactions. The parton-parton interaction cross section 
in the string-melting version of the AMPT model was 
taken to be 3 mb. 

The V2(pt) values obtained from all of these mod- 
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FIG. 21. (Color online) The difference in the vi values be- 
tween a particle X and its corresponding anti-particle (X) 
(see legend) as a function of ^/sjvjv a) and b) for 0-80% 
central Au+Au collisions. The dashed lines in plot a) are fits 
using Eq. (|18l) . and lines through the origin are shown for plot 
b). The values of /xb are from the parametrization of Ref. [4£| 
(see text for details). 



els were nearly identical for 7r + and 7r~. and K + and 
K~ , respectively. Only the anti-protons, compared to 
the protons, showed a significantly larger V2{pt) in the 
UrQMD and AMPT-SM models. This is in clear contra- 
diction to the observations from the data described here. 
The UrQMD model generally under-predicts the V2(pt) 
values. Only at 7.7 and 11.5 GeV are the anti-proton 
V2(pt) values close to or below the UrQMD values. As 



was pointed out above, a purely hadronic system (as de- 
scribed by the UrQMD model) does not appear to explain 
the relatively large flow of the particles at these energies. 
As seen in Fig. the AMPT-SM model provides the 
best description of the data, except for p and p at 7.7 
GeV. In all other cases, the AMPT default calculations, 
and more so the UrQMD calculations, under-predict the 
v2{pt) values. 



B. Interpretations from models 

Several interpretations have been suggested for the 
possible physical causes for the difference in the i>2 val- 
ues for particles and their corresponding anti-particles 
based on preliminary results. The process involved was 
to create or modify a model to qualitatively describe the 
difference in the V2 values between particles and cor- 
responding anti-particles that is shown in Fig. 1211 In 
Ref. [lH , it was argued that the effect results from quark 
transport from the projectile nucleons to mid-rapidity. 
The authors assumed that the elliptic flow of transported 
quarks is larger than that from produced quarks. Thus, 
the asymmetry of quarks and anti-quarks in the particles 
and corresponding anti-particles leads to a systematically 
larger flow of the particles compared to the anti-particles. 
The energy dependence was explained by the increase of 
nuclear stopping in heavy ion collisions with decreasing 
energy. The resulting patterns for 7T, K, p and A is qual- 
itatively in agreement with the data. However, a similar 
difference in v<i for mesons (n^, K^) and for baryons 

{{jp,p), (A, A), and (S~, is observed which is not 

expected in this picture. 

In Ref. [HI , an AMPT model calculation for ^/snn = 
7.7, 11.5 and 39 GeV was presented. The authors in- 
cluded mean-field potentials in the hadronic stage of that 
model. As a consequence of these potentials, particles 
like K~ and p are attracted by the hadronic matter and 
are trapped in the system whereas K + and protons feel 
a repulsive force and have the tendency to leave the sys- 
tem along the participant plane. The observed pattern 
shown in Fig. [2T1 cannot be explained by a default AMPT 
calculation without hadronic potentials, as discussed in 
Section [VII Al With the potentials included, a fair qual- 
itative agreement was achieved. However, the difference 
in i>2 between K + and K~ in the calculation is close to 
the difference for p and p, in clear contradiction to the 
present experimental results. The authors noted that 
further investigations are important to understand these 
effects in more detail. 

Similar studies were performed for the data collected 
by the KaoS collaboration at SIS (Schwerionensynchro- 
ton at GSI) at energies of 1-2 AGeV for K ± mesons [5(| . 
In this case, the IQMD (Isospin Quantum Molecular Dy- 
namics) transport model was used for the comparison to 
the KaoS results. The trends for the data and model 
calculations observed at those (very low) energies are op- 
posite to those reported in this paper. 
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FIG. 22. (Color online) The elliptic flow, 1)2, of ir , K , p and p as a function of the transverse momentum, pr, for 0-80% 
central Au+Au collisions at ^Jsnn = 7.7, 11.5 and 39 GeV. The symbols depict the data while the lines show the model results 
from AMPT with default settings (blue), AMPT with the string melting (SM) option and a hadronic cross section of 3 mb 
(red), and from UrQMD (black). The solid and dashed lines represent positively and negatively charged particles, respectively. 



In Ref. [52[, a hybrid (hydrodynamical plus UrQMD) 
calculation was performed. Qualitatively, the trend for 
At>2 can be described for protons, A, and 7r, whereas the 
trend for kaons is opposite to the present observations. 
The effect for the protons primarily resulted from the 
treatment of a non-zero net baryon number density and 
chemical potential. The results are slightly changed by 
using the UrQMD afterburner which describes the final 
stage interactions. Another effect discussed in this pa- 
per 52j is related to the event plane calculation. It was 
claimed that fluctuations in this calculation can bias the 
event plane to be rotated towards the most abundantly 
produced particles. This would, for example, increase the 
V2 values for protons and reduce them for anti-protons. 

In Fig. 1231 a study to explore this possibility is pre- 
sented. The elliptic flow for protons and anti-protons 
as a function of pr, for 0-80% central Au+Au collisions 
at y/sNN = 19.6 GeV is shown for two different kinds 
of reconstructed event planes. The event plane recon- 
structed using all of the charged particles is denoted by 
(+,-) EP, while the event plane reconstructed using only 
the negatively-charged particles is denoted by (-) EP. The 



v 2 (pt) values for protons using the (-) EP method are 
slightly, but systematically larger than those from the 
standard (+,-) EP method. The anti-proton V2{pt) val- 
ues are essentially unchanged. A reduced V2(pt) value 
would be expected for protons if such baryon number 
fluctuations caused such a bias. The increased V2{pt) for 
protons may be due to non-flow. For example, resonance 
decays could cause a larger change in non-flow contribu- 
tions to proton V2 than to anti-proton i> 2 between the two 
different event planes. More detailed studies from theory 
and experiment are needed to investigate the event-by- 
event baryon fluctuations and their possible effects on the 
event plane reconstruction. 

A recent calculation based on the Nambu- Jona-Lasinio 
(NJL) model can also qualitatively explain the differ- 
ences between p-p, A-A, and K + -K~ using the vector 
mean field potential [57j . The vector potential is repul- 
sive for quarks and attractive for anti-quarks, which re- 
sults in different flow patterns. To calculate the flow for 
the hadrons, a coalescence model was used. 

A different approach was followed in Ref. [58| by as- 
suming simplified rapidity distributions for u- and d- 
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FIG. 23. (Color online) The proton and anti-proton elliptic 
flow for 0-80% central Au+Au collisions at tJsnn = 19.6 
GeV, where "(+,-) EP" refers to the event plane reconstructed 
using all of the charged particles, and "(-) EP" refers to the 
event plane reconstructed using only the negatively charged 
particles. The error bars are statistical only. 



quarks that are different from those for s, u, d, and s 
quarks. It is claimed that under these initial conditions 
a breakdown of the v 2 NCQ scaling would not necessarily 
be a consequence of a phase transition, but rather the re- 
sult of the different rapidit y d istributions of the valence 
and produced quarks. Ref. [58[ also notes that the model 
results in a difference between particles and anti-particlcs 
that is opposite to that presented here. 



scaling with decreasing energy could be interpreted as the 
emerging dominance of hadronic interactions over par- 
tonic interactions in the systems formed in the collisions. 
The observed difference in the v 2 values demonstrates 
that the particles and anti-particles are no longer con- 
sistent with a single NCQ scaling law. The additional 
splitting between the particles and corresponding anti- 
particles at the lower beam energies breaks NCQ scal- 
ing. Even amongst the particles and anti-particles sep- 
arately, an absence of the baryon-mcson splitting is ob- 
served at y?SNN = 7.7 and 11.5 GeV in the representation 
vi(jtit — wo) for (% — mo) > 1 GeV/c 2 . However, the 
corresponding NCQ scaling shows no significant devia- 
tion from the scaling in the appropriate intermediate px 
range. 

It is observed that mesons at ^snn = 7.7 and 11.5 
GeV indicate a different trend at the highest pr values. 
This would be in agreement with the picture that the <j) 
mesons have a lower v 2 (pt) in a hadronic environment 
compared to other hadrons due to their lower hadronic 
cross section (37i - l40j . Larger event samples are needed 
at these energies in order to make more quantitative con- 
clusions. The corresponding anti-particles show a similar 
NCQ scaling trend as the particles at energies larger than 
ysKT^ = 11.5 GeV. At y/sWJJ = 7.7 and 11.5 GeV, the 
event sample sizes for most of the anti-baryons need to 
be increased in order to make quantitative statements on 
the validity of NCQ scaling. 

At energies larger than ^/snn = 11.5 GeV, NCQ scal- 
ing holds independently for particles and anti-particles, 
while at lower energies significant differences appear. The 
strong increase of the difference in v 2 between the parti- 
cles and corresponding anti-particles with decreasing en- 
ergy warrants further experimental and theoretical inves- 
tigation. 



C. Conclusions 

The strong energy dependence of the difference in 
v 2 {pt) between particles and their corresponding anti- 
particlcs is a new observation in the field of rclativistic 
heavy ion collisions. It cannot be explained in a purely 
hydrodynamic approach since particles and anti-particles 
have the same mass. It is also incompatible with a sce- 
nario where the flow is only produced in a thermalized 
and equilibrated QGP without any additional quark po- 
tentials. Other effects, such as hadronic interactions, 
or the scenarios involving transported quarks that were 
discussed in Section IVIIB1 could be responsible for the 
present observations. However, the agreement of the data 
with the transport based models is at present only quali- 
tative. The energy dependence of v 2 (X) — v 2 ( X) suggests 
a strong dependence on the values of the baryon chemical 
potential ps- 

The NCQ scaling was observed at ^fs^N = 200 GeV 
for all particles and anti-particles as they have the same 
number of quarks. The observed breakdown of such a 



VIII. SUMMARY 

Results on the mid-rapidity elliptic flow v 2 (px) for tt , 
K ± , K°, p, p, <f>, A, A, , S + , 0~ and Ti + from Au+Au 
collisions at y/sjyN = 7.7, 11.5, 19.6, 27, 39 and 62.4 
GeV were presented. For all of the particle species, v 2 
increases with increasing energy at high transverse mo- 
menta, whereas v 2 at low px values depend on the parti- 
cle species. A significant difference in v 2 (jpt) between 
the particles and corresponding anti-particles was ob- 
served. At energies above 39 GeV, the difference was 
approximately constant with energy, while the difference 
increased as the energy decreased. Hence, a significant 
dependence of v 2 {X) — v 2 (X) on the baryon chemical po- 
tential, [is, is indicated. The difference v 2 (X) — v 2 (X) 
was larger for baryons than for mesons. This difference 
cannot be reproduced by transport models in their stan- 
dard configuration. Only the AMPT model with an in- 
cluded hadronic potential shows a similar pattern. Other 
models show the same qualitative trends. 

The NCQ scaling that was observed for all particles 



and anti-particles at ^snn = 200 GcV, no longer holds 
at the lower beam energies of ^snn = H-5 and 7.7 GeV. 
This is seen as an increase of V2{X) — i>2(X) with de- 
creasing beam energy. The baryon-meson splitting of 
v-iirriT — rrio) for (% — mo) > 1 GeV/c 2 , which formed 
the basis of NCQ scaling observation at 200 GeV, was not 
observed for anti-particles at the lower energies. In the 
representation of V2 (mx — mo)/n q vs. (jut — mo)/n q , no 
significant deviations from NCQ scaling were observed for 
particles and anti-particles separately at energies above 
^/sWW = 11.5 GeV. At y/smf = 7.7 and 11.5 GcV, the 
anti-baryons and the </> meson indicate a different trend. 
At the highest (mr-mo) data points at -^/ijvw = 7.7 and 
11.5 GeV, these particles deviate from the other hadrons 
by 1.8ct and 2.3<t, respectively. 
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APPENDIX: NCQ SCALING OF v 2 (p T ) 



The NCQ scaling of v 2 is shown in Fig. [24] Plotted there is v-2/n q versus the scaled transverse momentum, pr/n q , 
where n q is the number of quarks in the particle. The same data is shown versus (mj- — m )/n q in Fig. 1191 There is 
a wider variation of the scaled v 2 when plotted versus PT/n q as compared to (m^ — m )/n q . 
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FIG. 24. (Color online) The number-of- Constituent Quark (NCQ) scaled elliptic flow, vifn q versus pr/uq, for 0-80% central 
Au+Au collisions for selected particles a) and corresponding anti-particles b). The dashed lines show the results of simultaneous 
fits with Eq. (|17(l to all particles except the pions. 

Simultaneous fits to all of the Vijn q versus pT/n q values, except for those for the pions, were performed. In 
similarity to Fig. [201 shown in Fig. (5S]is the ratio of the data points in Fig. [M]to the simultaneous fits as a function 
oipT/n q . Most of the data points are within 10% of the fit function at pr/nq values larger than 1 GeV/c. At lower 
momenta, the V2/n q versus pr/nq values diverge due to the mass splitting that was shown in Fig. [TOl 
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FIG. 25. (Color online) The Number-of-Constituent Quark (NCQ) scaled elliptic flow, V2(pT)/n q , ratio to a fit function (see 
text) for 0-80% central Au+Au collisions for selected particles a) and corresponding anti-particles b). Most of the data points 
at the larger pr/riq values are consistent with unity to ±10%, which is shown as the shaded areas to guide the eye. Some of 
the data points for <f> and S are outside of the plot range. 



33 



REFERENCES 



M. M. Aggarwal et 



[STAR Collaboration], 
Phys. Lett. B 507, 121 



arXiv: 1007.2613 [nucl-ex]. 
D. Kharzeev and M. Nardi, 
(2001). 

R. Snellings, New J. Phys. 13, 055008 (2011). 

S. Voloshin and Y. Zhang, Z. Phys. C 70, 665 (1996). 

A. M. Poskanzer and S. A. Voloshin, Phys. Rev. C 58, 
1671 (1998). 

J. -Y. Ollitrault, Phys. Rev. D 46, 229 (1992). 

H. Sorge, Phys. Rev. Lett. 82, 2048 (1999). 

S. A. Voloshin, A. M. Poskanzer and R. Snellings, in 

Landolt-Boernstein, Relativistic Heavy Ion Physics, Vol. 

1/23, p. 5-54 (Springer- Verlag, 2010). larXiv:0809.2949l 

[nucl-ex] . 

N. Borghini and J. Y. Ollitrault, Phys. Rev. C 70, 064905 

(2004). 

P. Sorensen. larXiv:0905. 01741 [nucl-ex] . 

J. Adams et al. [STAR Collaboration], Nucl. Phys. A 

757, 102 (2005). 

K. Adcox et al. [PHENIX Collaboration], Nucl. Phys. A 
757, 184 (2005). 

P. Huovinen, P. F. Kolb, U. W. Heinz, P. V. Ruuskanen 
and S. A. Voloshin, Phys. Lett. B 503, 58 (2001). 
C. Nonaka, R. J. Fries and S. A. Bass, Phys. Lett. B 583, 
73 (2004). 

T. Hirano and Y. Nara, Phys. Rev. C 69, 034908 (2004). 
C. Shen and U. Heinz, Phys. Rev. C 85, 054902 (2012) 
[Erratum ibid. C 86, 049903 (2012)] 

B. Schenke, S. Jeon and C. Gale, Phys. Rev. Lett. 106, 
042301 (2011). 

C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 
87, 182301 (2001). 

J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 
92, 052302 (2004). 

J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 
95, 122301 (2005). 

B. I. Abelev et al. [STAR Collaboration], Phys. Rev. C 
75, 054906 (2007). 

B. I. Abelev et al. [STAR Collaboration], Phys. Rev. C 
81, 044902 (2010). 

B. I. Abelev et al. [STAR Collaboration], Phys. Rev. Lett. 
99, 112301 (2007). 

B. Mohanty and N. Xu, J. Phys. G 36, 064022 (2009); K. 
J. Wu, F. Liu and N. Xu, J. Phys. G 37, 094029 (2010). 
L. Adamczyk et al. [STAR Collaboration], Phys. Rev. C 
86, 054908 (2012). 

K. H. Ackermann et al., Nucl. Instr. and Meth. A 499, 
624 (2003). 

M. Anderson et al, Nucl. Instrum. Meth. A 499, 659 
(2003). 

W. J. Llope (STAR TOF Group), Nucl. Instr. and Meth. 
B 241, 306 (2005). 



[29] H. Bichsel, Nucl. Instr. and Meth. A 5 62, 154 (2006). 
[30] H. Masui and A. Schmah, larXiv:1212.3650l [nucl-ex] 
(2012). 

[31] B. I. Abelev et al. [STAR Collaboration], Phys. Rev. C 

77, 054901 (2008). 
[32] F. Wang, M. Nahrgang and M. Bleicher, Phys. Rev. C 

85, 031902 (2012). 
[33] X. Dong, S. Esumi, P. Sorensen, N. Xu and Z. Xu, Phys. 

Lett. B 597, 328 (2004). 
[34] P. Huovinen and P. V. Ruuskanen, Ann. Rev. Nucl. Part. 

Sci. 56, 163 (2006). 
[35] J. Beringer et al. [Particle Data Group Collaboration], 

Phys. Rev. D 86, 010001 (2012). 
[36] Y. Burnier, D. E. Kharzeev, J. Liao and H. -U. Yee, Phys. 

Rev. Lett. 107, 052303 (2011). 
[37] A. Shor, Phys. Rev. Lett. 54, 1122 (1985). 
[38] A. Sibirtsev, H. -W. Hammer, U. -G. Meissner and 

A. W. Thomas, Eur. Phys. J. A 29, 209 (2006). 
[39] H. van Hecke, H. Sorge and N. Xu, Phys. Rev. Lett. 81, 

5764 (1998). 

[40] Y. Cheng, F. Liu, Z. Liu, K. Schweda and N. Xu, Phys. 

Rev. C 68, 034910 (2003). 
[41] D. Adamova et al. [CERES Collaboration], Nucl. Phys. 

A 894, 41 (2012). 
[42] S. A. Voloshin, Nucl. Phys. A 715, 379 (2003). 
[43] A. Adare et al. [PHENIX Collaboration], Phys. Rev. C 

85, 064914 (2012). 

[44] R. J. Fries, B. Muller, C. Nonaka and S. A. Bass, Phys. 

Rev. C 68, 044902 (2003). 
[45] B. Mohanty and N. Xu, J. Phys. G 36, 064022 (2009). 
[46] M. .Nasim, B. Mohanty and N. Xu, larXiv:1301.1375l 

[nucl-ex] . 

[47] L. Adamczyk et al. [STAR Collaboration], 

larXiv: 1212.33041 [nucl-ex] . 
[48] S. K. Tiwari, P. K. Srivastava and C. P. Singh, Phys. 

Rev. C 85, 014908 (2012). 
[49] B. I. Abelev et al. [STAR Collaboration], Phys. Rev. C 

79, 034909 (2009). 
[50] S. A. Bass et al, Prog. Part. Nucl. Phys. 41, 255 (1998); 

M. Bleicher et al, J. Phys. G 25, 1859 (1999). 
[51] Z.-W. Lin et al, Phys. Rev. C 72, 064901 (2005). 
[52] J. Steinheimer, V. Koch and M. Bleicher, Phys. Rev. C 

86, 044903 (2012). 

[53] Wang X N and Gyulassy M 1991 Phys. Rev. D 44, 3501 
[54] J. C. Dunlop, M. A. Lisa and P. Sorensen, Phys. Rev. C 

84, 044914 (2011). 
[55] J. Xu, L. -W. Chen, C. M. Ko and Z. -W. Lin, Phys. 

Rev. C 85, 041901 (2012). 
[56] F. Uhlig et al., Phys. Rev. Lett. 95, 012301 (2005). 
[57] T. Song, S. Plu mari, V. Greco, C. M. Ko and F. Li, 

larXiv:1211.55TT1 [nucl-th] . 

[58] V. Greco, M. Mitrovski and G. Torrieri. larXiv: 1201.48001 

[nucl-th] . 



